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Abstract 

We investigate the distribution of oscillator strengths for the recombination of excitons 
in a two dimensional sample, trapped in local minima of the confinement potential: the 
results are derived from a statistical topographic model of the potential. The predicted 
distribution of oscillator strengths is very different from the Porter-Thomas distribution 
which usually characterises disordered systems, and is notable for the fact that small 
oscillator strengths are extremely rare. 



1. Introduction 

Many experiments have been performed in which excitons (excitations in the form 
of an electron-hole bound state) are trapped in a layer of lower bandgap semiconductor, 
between two layers of higher bandgap material [1,2]. Exciton absorption and emission lines 
in these two-dimensional samples are typically much broader than in three-dimensional 
systems. This is often due to fluctuations in the width of the layer: there is a 'quantum 
confinement' contribution to the energy of the exciton, which is a decreasing function E{w) 
of the layer width w, and which is analogous to the ground state energy tt'^H^ /2mw'^ of a 
particle trapped in a one-dimensional potential well [3,4]. 

The excitons are able to interact with phonons by radiationless processes [5,6], be- 
coming trapped in local minima of the confinement energy: the lifetime for electron- hole 
recombination is typically long enough that most of the luminescence observed from these 
samples is from trapped excitons. This trapping effect has recently been observed directly 
by spatially resolved studies of exciton luminescence [7,8]: if the luminescence is recorded 
from a macroscopic area of the sample, a broad spectrum is observed, whereas sufficiently 
small microscopic areas show either no luminescence, or a small number of relatively sharp 
lines, corresponding to excitons in a single trap. These exciton traps may be thought of a 
novel type of mesoscopic system. 

The spectra of two dimensional excitons show some near-universal features, for ex- 
ample it was noticed that the Stokes shift S of the luminescence spectrum relative to the 
absorption peak, and the width W of the absorption spectrum usually satisfy S/W ~ 0.6, 
independent of the semiconductor materials or of the magnitude of the broadening [9]. 
This observation was explained using a classical picture of trapping of excitons in local 
minima of a smooth effective confinement potential, modelled as a Gaussian random func- 
tion [10]. This statistical topographic model predicted S/W — 0.55.., in good agreement 
with the experimental values. In section 2 we review this model, and discuss a variant of 
the experimental approach which may give better agreement with theoretical predictions. 

In this paper we discuss the distribution of oscillator strengths (or, equivalently, Ein- 
stein A coefficients) for exciton recombination, for excitons localised in minima of the 
effective confinement potential. We use the same statistical topographic model as was 
previously used to calculate the Stokes shift. If the model is correct, the probability dis- 
tribution of oscillator strengths for localised excitons should be a near universal signature 
of trapped excitons, but we believe that this may be a more sensitive test of the validity 
of the model. The predicted probability distribution of oscillator strengths is 

P{I) = exp(l/2/^)erfc(V3/2/2) (1.1) 

(the scaling of the intensities / is arbitrary; for this form of the distribution the mean 
intensity is (/) ~ 0.699..). This result is interesting because the predicted distribution has 
very few small intensities: the fraction of recombination lines predicted to be less than 
half the mean intensity is approximately 3.15 x 10~^. It is very different from the Porter- 
Thomas distribution, P{I) ~ exp(— /)/\/7, which usually characterises the distribution 
of oscillator strengths for disordered or complex systems [11], and which follows from the 
Gaussian distribution of matrix elements in such systems. 

Section 3 discusses the model used for calculating the matrix elements for recombina- 
tion of excitons trapped in local mimima of a potential. In section 4 we give a complete 
characterisation of the distribution of quadratic forms characterising the stationary points 
of an isotropic Gaussian random function: although this is not difficult to obtain using 
standard techniques of statistical topography [12,13], we could not find a discussion of it 
in the literature in the form we require. These results are then used to obtain (1.1). 



The field of statistical topography was largely stimulated by the desire to understand 
optical properties of random surfaces, such as that of the sea. In section 4 we also comment 
on the closest optical analogy of our results, namely the distribution of intensities of 
reflected images of a small light source on a distant random surface (for example, 'sea 
glitter', reflections of the sun from the sea, observed from an aircraft [14]). The distribution 
(1.1) turns out to be very different from that of the intensities of sea glitter sparkles. 

2. The statistical topographic model 

Here we briefly review the statistical topographic model for exciton luminescence 
presented in [9,10]. We also propose a reason for the small discrepancy between theory 
and experiment, and a variant of the experimental approach which may give a better 
agreement with experiment. 

Figure 1 shows the absorption (a) and luminescence (b) spectra of excitons in a semi- 
conductor heterostructure: the data are taken from [15]. First consider the form of the 
absorption peak. In the absence of disorder the peak would be very sharp, because (un- 
like unbound electron-hole pairs) conservation of momentum implies that the oscillator 
strength for creation of the exciton vanishes unless the centre of mass of the exciton is 
stationary ([16]; this can also be verifled using the model discussed in section 3). The 
width of the absorption peak in flgure 1(a) is determined by inhomogeneous broadening 
due to disorder. The energy of an exciton in the ground state depends upon the width of 
the potential well in which it is confined, and the well width varies randomly with position 
in the plane. We assume that the length scale over which the well width varies is large 
compared to the width of the wells: this implies that the ground state energy of a static 
exciton at position {x, y) is a well defined smooth function, which we denote by E{x^ y). It 
is reasonable to assume that the fiuctuations of E[x, y) represent contributions from many 
independent events, and the central limit concept then indicates that the fiuctuations are 
Gaussian. Because the excitons are created with zero centre of mass motion, the absorp- 
tion spectrum is proportional to the distribution of E{x, y). The curve (c) in figure 1 is a 
fit of a Gaussian curve to the absorption spectrum (a): it fits quite closely. 

Now consider the luminescence peak. This is shifted toward lower energies because 
the exitons can lose energy before they decay. Time-resolved studies of spectral hole 
burning indicate that the energies of excitons can change over a timescale of typically 
tens of picoseconds, much shorter than the half-life for decay of excitons, typically several 
hundred picoseconds [5,6]. The predominant mechanism of energy loss for the excitons 
appears to be by the excitation of phonons: if the absorption spectrum is probed with 
narrow spectral lines, it is possible to observe features in the luminescence spectrum which 
are shifted from the probe frequency by multiples of the frequency of the optical phonons 
[6]. These results justify the following picture of the luminescence process: after the exciton 
is created at position (xo^yo) with energy Eq — E{xo,yo), it will move into regions where 
the potential energy E{x, y) is less than E'o, and the excess energy Eq — E{x, y) appears 
as kinetic energy. The moving exciton is able to excite phonons, and as it does so it loses 
kinetic energy. Eventually, if it does not decay in the meantime, it will end up trapped 
in a local minimum of the potential energy E{x, y). Because the exciton lifetime is much 
longer than the timescale associated with energy transfer to phonons, most of the excitons 
are trapped close to a local minimum of the potential energy E{x^ y) when they decay. 

These considerations lead to a model in which both the absorption and luminescence 
spectra are determined by the statistical topography of a Gauss random function, which 
we will denote by f{x,y), and which represents the energy function E{x, y) after applying 



linear scaling transformations to E, x and y such that 



{fl) = ify) = 1 • (2.1) 

The absorption spectrum is proportional to the Gaussian probability distribution of this 
function, whereas the luminescence spectrum is determined by the distribution of heights 
of its local minima. 

The Gaussian random function f{x,y) is characterised completely by its correlation 
function, C(r), which in the isotropic case is a function of r = |r| only. The calculations 
in [13], [10] show that the distribution of stationary points depends on the correlation 
function only through averages of second derivatives. In the isotropic case, after scaling 
the function so that (2.1) is satisfied, the distribution of stationary points is characterised 
by a single parameter a [10]: 

ifL) = ifyy) = HO = 3a . (2.2) 

The distribution of heights of local minima was determined analytically in [9,10] : the result 
is 
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where erfc(a;) is the complementary error function [17]. 

The prediction (2.3) contains the imdetermined parameter a. In [9,10] it was ar- 
gued that the annealing process involved in the growth of the heterostructures causes the 
fluctuations of E{x,y) to be suppressed by a diffusive process: we write 

E{x,y)^ J dx' J dy'P{x-x',y-y')E'{x',y') 

= 8^ exp[-(X2 + Y')/4Dt] (2.4) 

where E'(x,y) characterises an initial distribution of well width fluctuations with a much 
shorter correlation length, D is the diffusion constant and t the annealing time. We may 
therefore model E[x^ y) as the convolution of a white noise function with a Gaussian. It 
follows that the correlation function of E[x, y) is also a Gaussian, implying that a = 1. 

The distribution of minima (2.3) with a = 1 is shown as curve (d) in figure 1: the mean 
and variance have been scaled to correspond to the mean and variance of the Gaussian 
distribution (c). The curve is not a particularly good fit to the luminescence spectrum, 
(b). The experimental data shown in figure 1 are typical: it is usually found that the 
Stokes shift is somewhat higher than our theoretical prediction [9]. We will make two 
points about this discrepancy. 

First, wc propose a qualitative explanation of this observation: the excitons can move 
into deeper local minima than the ones in which they were initially trapped before they 



recombine, either by thermally assisted hopping or quantum mechanical tunnelling. This 
explanation is supported by the literature on time resolved luminescence studies, which 
show that the Stokes shift initially assumes a value close to our prediction, and then slowly 
increases [18]. 

Our second, more important, point concerns resonant Rayleigh scattering experiments 
on excitons, such as that discussed in [15]. /,From the discussion above, it is apparent that 
it would be desirable to have a more direct probe of the density of local minima of the 
effective potential, to facilitate comparison between theory and experiment. We will now 
argue that resonant Rayleigh scattering spectra measure the density of local minima of the 
effective potential. The intensity of Rayleigh scattering from an exciton will be proportional 
to the time over which the exciton survives at the energy of the incident radiation. The 
lifetime of an exciton trapped in a minimum of the effective potential is equal to the 
lifetime for luminescent decay, typically less than a nanosecond. An exciton which is not 
trapped may lose energy by exciting phonons, and the lifetime for these processes is much 
shorter, typically a few tens of picoseconds. The spectra for resonant Rayleigh scattering 
are proportional to a density of states weighted by the exciton survival time. Because the 
survival time is very much larger for trapped excitons, the Rayleigh scattering spectrum 
is weighted very heavily by the density of states for trapped excitons. 

The isolated points plotted on figure 1 are resonant Rayleigh scattering amplitudes 
measured at a discrete set of frequencies, taken from data published in [15]. Both the mean 
and the variance of these data are remarkably close to the theoretical distribution of local 
minima, curve (d), which is a plot of (2.3) with a = 1, and with the mean and variance 
chosen to match those of the Gaussian fit to the absorption spectrum. 

3. Model for exciton luminescence matrix elements 

We consider a solid consisting of N atoms, each of which can exist in four different 
states: 

1. Neutral and unexcited, |0). 

2. With additional electron (in the conduction band), |— ). 

3. With an electron removed (i.e. with hole in valence band), |+). 

4. Excited atom (or Frenkel exciton), with an electron in the conduction band and a hole 
in the valence band, jci?) = |±). 

We consider a Hilbert space for the solid with 4^ basis vectors, consisting of all possible 
combinations of these four states of the individual atoms: for example a state with two 
electrons, a hole and an exciton could be written |0, 0, — , 0, 0, .., 0, 0, +, 0, ±, 0, — , 0, ..) = 
|ni, n2, ns; ns, n^), where ni, n2 are the positions of the electrons, is the position of the 
hole, and ns is the position of the Prenkel exciton. 

This model is reasonable for insulators in which the electrons are tightly bound to 
individual atoms. This description can also be carried over directly to semiconductors, if 
the electron and hole states localised on individual atoms are replaced by Wannier states 
derived by integrating over the conduction and valence band wavefunctions respectively 



In a semiconductor, the excitons are typically of the Wannier type, in which the 
electrons and holes, although correlated, are not bound to the same atomic orbital. The 
Wannier exciton state \ew) is a superposition of electron and hole states of the form 



[16]. 




(3.1) 
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where rii, n2 are the postions of the electron and hole respectively. In the case of the weakly 

localised Wannier cxciton, the coefficients c^^.na can be approximated by a continuous 
wavefunction ?/;(x, y), where x and y are the locations of atoms with labels ni and n2. 
The wavefunction ip{:x., y) is an eigenfunction of an effective Hamiltonian [16] 

H = :^pI + :^pI + , ^ , + V,{r,) + ]/h(rh) • (3.2) 
2me 2mh 47reo|re-rh| 

The intensity of emission from an exciton state is proportional to the square of the 
dipole matrix element (cvkI-^IO) for the transition between the exciton state and the ground 

state of the system, summed over three orthogonal choices for the coordinate X. This 
matrix element can be expressed in terms of the dipole matrix elements of the localised 
basis states. We assume that the dipole matrix element for the basis states is 

{0\X\ni;n2) = e 6nin2 (3-3) 



i.e. the dipole matrix element for recombination is assumed to be negligible, unless the 
electron and hole are on the same atomic site. Most crystals have centres of symmetry, and 
the dipole matrix element (0|X|n;n) vanishes if the atomic orbitals or Wannier functions 
defining the valence and conduction bands have the same parity with respect to the X 
coordinate. The atomic orbitals associated with successive bands typically have opposite 
parity with respect to one of the coordinates, so that for at least one choice of the co- 
ordinate X the matrix element considered in (3.3) does not vanish because of symmetry 
considerations. 

The required matrix element can now be calculated using (3.3): 

(0|X|ni;n2) 

nin2 

= e'^Cnn^ I dr V'(r, r) (3.4) 

n 

(here e' is another constant). The matrix element is therefore proportional to the amplitude 
for the electron and hole to be at the same site. 

We now consider the case of excitons in heterostructures, where the excitons are 
trapped in a layer of low bandgap material between regions of higher bandgap. Imper- 
fections of the growth process result in random fluctuations of the layer width: we will 
assume that the fluctuations of the layer width are on a larger scale than the size of the 
exciton [10], and we therefore use the following model 

V;(re) + Vi(rh) = ^^e(^e) + V^{Zh) + We(Xe) + Wh(Xh) 

~ ^'e(^e) +t'h(^h) + We(X) -l-tth(X) (3.5) 

where Xg/^ = (iCe/hjye/h) cire the electron/hole coordinates in the (x, y) plane, and X, 
x are the corresponding centre of mass and relative coordinates. In the second line the 
assumption that the variation of the potentials -Ug/h is slow on the scale of the exciton 
diameter justifles the approximation Xg ~ Xh ~ X. The effective Hamiltonian can then be 
written in the separated form 



(re, r/i, Pe, Ph) = Hq^-Hi 



1 

Ho = — + r-r + Ve{Ze) + Vh{Zh) 

2// 47reo|r| 
^i = ^P' + KfT(X) (3.6) 

where Veff(X, Y) = Ue{X, Y) + 'Uh(X, F), P is the momentum conjugate to X = (X, Y), 
and |U, M are the reduced and total masses. The corresponding solution of the Schrodinger 
equation is a product of an exciton wavefunction xi.^^ V-i ^e-, ^h), and a wavefunction 4){X, Y) 
for the centre of mass motion in the (X, 1") plane. Equations (3.4) and (3.6) show that 
the transition strength for exciton recombination is of the form 



I ^ \{0\X\ew)f = C 



dX J dY (f>{X, Y) 



(3.7) 



where the constant C is the same for all exciton states. 

We first discuss the interpretation of this result for a non-disordered system. In this 
case, the centre of mass wavefunctions are (/>(X) = exp[ik.X], where A is the area 

of the sample. The recombination transitions only occur from the ground state, for which 
(f){X, Y) = A~^/'^, because the integral in (3.7) vanishes for all of the other states. Equation 
(3.7) then implies that the transition rate for recombination from the centre of mass ground 
state is proportional to the area of the system, whereas the rate for all of the other possible 
states is zero. This is in accord with the expectation that the sum of the transition rates 
for a system of area A should be proportional to A. 

Now consider the case of a disordered system. We will use the same model as in [9] and 
[10]: we assume that the excitons interact strongly with phonons, and that mobile excitons 
rapidly lose energy by exciting phonons: most of the exciton recombination therefore occurs 
after the excitons have become trapped in minima of the effective potential Ves{X, Y). We 
will assume that the effective potential for the centre of mass motion is an isotropic Gauss 
random function. The form of the effective potential in the neighbourhood of its minima 
can be approximated by a quadratic form: if the origin of the {X, Y) plane is shifted to lie 
at the minimum, we write 

Fefr(^, Y)^Vo + ^[Vxx X^ + Vyy + 2Vxy XY] ^Vo + ^X'^MX (3.8) 

where M is the Hessian matrix of second derivatives evaluated at the minimum. The 
exciton recombination occurs from a ground state of the centre of mass motion trapped in 
this minimum, for which the wavefunction (j){X, Y) is a harmonic oscillator ground state, 
which satisfies 

J dX J dY (p{X,Y) = c[detM]-^/^ (3.9) 

where c is independent of the local environment in which the exciton is trapped. The 
distribution of transition strengths / for trapped excitons is therefore determined from the 
distribution of determinants of the Hessian matrix at minima: 



Ir^\D\-'^/^, D = det{M) 



(3.10) 



4. Distributions of intensities 



4.1 Distribution of quadratic forms at statioary points 

Methods for calculating properties of point singularities such as minima are well 
known: the one and two dimensional cases are discussed by Rice [12] and Longuet-Higgins 
[13]. We will calculate the distribution of the trace T and determinant D of the Hessian 
matrix 

fxx fxy 



M = 



fxy fyy 



(4.1) 



describing the second derivatives fxx, fxy, fyy of an isotropic Gauss random function f{x,y) 
at its stationary points. 

We follow the approach and notation of [10]: the distribution of extrema is determined 
by the joint probability distribution P{f, fx, fy, fxx, fyy, fxy) of the function /, its first 
derivatives fx, fy, and its second derivatives evaluated at the same point [x, y). By a simple 
adaptation of the calculation in [10], the joint distibution of the trace and determinant is 

P(T,D) = -1 Jdf J dfxxj dfyyj df yyP{f, 0,0, fxx, fyy, fxy)\d\5{D-d)S{T-t) (4.2) 

where is a normalisation factor, and t = {fxx + fyy), d — fxxfyy — fxy The function is 
assumed to be scaled so that (3.1) is satisfied. 

In terms of the variables R = ^t, X = ^{fxx — fyy), Y = fxy, the probability 
distribution in (4.2) was obtained in [10]: 



P{f, 0; 0; fxx, fyy, fxy) — 



:exp(-/V2) 



(27r)3aV2a - 1 

X exp[-(/ + i2)V2(2a - 1)] exp[-{X^ + Y^)/2a\ . (4.3) 
The integral (4.2) can now be evaluated easily: the normalised probability density is 



P{T,D) 



1 



2,2a? 




\D\ exp[i:>/2a] e-ycp[-ZT'^ /lQa]Q{D -T'^ /A) 



(4.4) 



where Q{x) is a step function, decreasing from 1 to at a; = 0. 

4.2 Distribution of oscillator strengths 

Now we will use the results of section 3 to calculate the distribution of transition 
strengths for excitons. The distribution of the determinant D of the Hessian matrix at 
stationary points is obtained by integrating (4.4) over T: 



P{D)= I dT P{T,D) = -^\D\ exp{D/2a) 



e{D) + Q{-D) erfc(|V3^/a) 



(4.5) 



It is noteworthy that the parameter a only appears in the ratio D/a in this probability 
measure: it therefore sets the scale of the distribution of D but does not alter its functional 
form. The distribution of intensities is therefore insensitive to this parameter, and we may 
take a = 1. 

The minima are those extrema for which D > and T > 0. The distribution of 
intensities / ~ follows immediately from this expression, considering only the 



branch with D > 0. The result, with arbitrary scale for the intensities, is given by (1.1). 
Numerically, the mean value of the intesities with distribution (1.1) is found to be (/) = 
0.699... The distribution (1.1) is plotted in figure 2(a). 

4.3 Relation to sea glitter 

The analysis of Gauss random functions in two dimensions was largely motivated by 
the desire to achieve a statistical understanding of the surface of the sea. It is natural to 
ask which experimentally observable property of the sea surface corresponds most closely 
to our distribution of exciton transition intensities, and to compare the two results. The 
corresponding property is 'sea glitter' [14] , which is the pattern of reflections of the sun on 
the ocean surface seen from a high flying aircraft. 

To simplify the dicussion we will assume that the sun is overhead, and that the sea 
glitter is observed by looking vertically downwards. If the radius of the sun subtends an 
angle e, the observer sees a bright area on the sea suface due to reflected sunlight whenever 
the angle of the sea surface is smaller than |e. There is therefore a bright spot on the sea 
surface in the neighbourhood of every maximinu, minimum or saddle point: each of these 
bright spots is seen as an ellipse, of area A ~ |det(M)|~^, where M is the Hessian matrix 
at the stationary point. The quantity which corresponds to the distribution of exciton 
intensities is the distribution of integrated intesities / ~ A of the reflected images of the 
sun. The distribution of determinants, for all types of extrema, is given by (4.5) and the 
corresponding distribution of / ~ \D\~'^ is 

P{I) = [exp(-l/2/) + exp(l/2/) eI^c{l^/3/I)] . (4.6) 

Note that, once again, the distribution is independent of the parameter a: this result is 
therefore universal for isotropic sea surfaces. This distribution is plotted in figure 2(b) 
(again, the normalisation of the intensity distribution is arbitrary, and in this case we find 
numerically (/) = 0.433..). This result is very different from the distribution of exciton 
intensities. The reasons for the difference are that all of the extrema contribute, and that 
there is a different relationship between / and D = det(M). 
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Figure Captions 



Figure 1. Experimental data from [15]: a) absorption spectrum, b) luminescence spectrum, 
isolated points - resonant Rayleigh scattering. Theoretical curves: c) Gaussian fit to 
absorption spectrum, d) corresponding distribution of heights of local minima. 

Figure 2. Predicted distribution of intensities of: a) exciton luminescence, b) sea glitter 
sparkles. 
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